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ABSTRACT: We have observed that the separation of variables method and graph theory in 


mathematics can be applied to several types of problems on bounded three-dimensional domains. 
The problem must be linear and have homogeneous boundary-conditions. In this paper, we present 
the result of classical problem in three dimensions by applying modified integral form, of Newton's 
law of cooling for sphere, which is encountered in ordinary differential equations texts. Also show 
that, the rate at which a body cools is proportional to the difference of its temperature and the 
temperature of the environment. The assumed symmetries in the problem will permit us to reduce 
the dimension of the problem to spatial dimension and time with the help of Bipartite Graphs in 
Graph theory. 
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INTRODUCTION 


The below graph is Bipartite Graphs - A 
simple graph is said to be bipartite if its vertex 
set can be divided into two disjoint sets such 
that every edge in has its initial vertex in the 
first set and the terminal vertex in the second 
set. Total number of edges are (n*m) with 
(n+m) vertices in bipartite graph. 


A simple graph is bipartite if and only if it is 
possible to assign one of two different colors to 
each vertex of the graph so that no two adjacent 
are assigned the same color. A bipartite graph 
with and vertices in its two disjoint subsets is 
said to be complete if there is an edge from 
every vertex in the first set to every vertex in 
the second set, for a total of edges[1]. 


The Meijer’s [2] 


G-function of one variable was defined in terms 
of Mellin-Barnes type integrals as follows: 
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i ane 
and an empty product is interpreted as 1,0 <m 
<q, 0 <n <p, and the parameters are such that 


no poles of I'(b; - s) (j = 1, ..., m) coincides with 
any pole of [(1 - aj +s) (j= 1, ..., n). 


There are three different paths ‘L’ of 
integration 


(i) ‘L’ runs from — i> to + ico so that all poles 
of I'(b; - s) (j = 1, .... m) are to the right 
and all the poles of ['(1 - aj +s) § =1,... 
n) to the left of L. The integral converges if 
p+q<2 (m+n) and jargx|<(m+n-¥% 
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-¥,q)r. 

(ii) ‘L’ is aloop starting and ending at + and 
encircling all poles of I'(b; - s) (j = 1, ..., m) 
once in the negative direction, but none of 
the poles of I(1 - aj + s) (j = 1, ...,n). The 
integral converges if q = 1 and either p <q 
or p=qand |x| <1. 

(iii) ‘L’ is a loop starting and ending at — © and 
encircling all poles of [(1 -aj+s) (j=1 
n) once in the positive direction, but none 
of the poles of I'(b; - s) (j = 1, ..., m). The 
integral converges if p => 1 and either p > q 
or p=qand |x| < 1. 

Given a sphere, whose initial temperature be 
contingent only on the distance from the 
epicenter (for instance- a constant initial 
temperature)[4] and whose boundary is kept at 
a persistent temperature; predict the 
temperature at any fact privileged the sphere at 
a later time. This is the problem, of determining 
the temperature of the center of a potato that 
has been put in a hot oven. The reader might 
conjecture that the given problem is imperative 
for medical field who, want to determine the 
time of an individual's death. Early researchers 
used this problem to regulate the age of the 
earth based on assumptions about its initial 
temperature and its final temperature [2]. 


This cooling-problem may remind the reader 
of Newton's law of cooling, which is 
encountered in ordinary differential equations 
texts. Recall, that this regulation status that the 
rate at which a body cools is proportional to 
the difference of its temperature and the 
temperature of the environment. Statistically, if 
T = T(t) is the temperature of a body and T, is 
the temperature of its environment, 


Then T'(t) = K(T--T), 
Where, K is the constant of proportionality. 


But the reader should point-out that this law 
applies only in the circumstance that the body 
has a uniform, homogeneous temperature. In 
the PDE queries are considering, the 
temperature may vary radically throughout the 
body. In this paper, we shall make application 
of following modified form of the integral [5] 
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J, (sin x)? sin nxdx = a EH (2) 
where, Re (s) > Oand Legendre’s duplication 
formula 
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Va (2z) = 272-41 (z) T(z + ¥). (2) 


2. INTEGRAL 


The integral, which we need as follows: 
Ty. wWw—l1a: Ge (a aj, 3 ie Pld 
f (sin x) sinnx Gyq |Z (sinx)? li, »,, x 


= m sin 


eqns [ tE- 
ia cans) 


Dro] c4 
1) | 


provided that m > 0, Re[w + 2min,<j<m bj] > 0, 


larg z| <(m+/-%p-% q)x[6] [7]. 


Proof: 


The integral (4) can be obtained easily by 
making use of the definition of G-function as 
given in (1) and the formulae (2) and (3). 


3. FORMULATE THE PROBLEM 


For simplicity, let us consider a sphere of 
radius p =7 whose initial temperature is To = 
constant. We will assume that the boundary is 
held at zero degrees for all time t >0. If u is the 
temperature, then in general u will depend on 
three spatial coordinates and time. But a little 
reflection shows that the temperature will 
depend only on the distance from the center of 
the sphere and on time [8]. Evidently, the 
temperature u must satisfy the heat equation 


ur= kAu, (5) 


where k is the diffusivity and A is the Laplacian. 
It should be clear that a spherical coordinate 
system p,@,@is more appropriate than a 
rectangular system, and u = u(p, t). Because u 
does not depend on the angles gand9Q, the 
Laplacian takes on a particularly simple form: 


Therefore, we may formulate the model as 


ut = k(uga + - up), 0 <p<z,t> 0. 
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u(p, 0) = To, O S p<z. 
(7) 


Observe that there is an implied boundary 
condition at p= 0, namely that the temperature 
should remain bounded. 


To solve problem (5)-(7) we assume u(p, t) = 
y(p)g(}[9]. Substituting into the PDE and 
separating variables gives, following solution, 
which is given in [9, p. 144-146]: 
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ulp,t) = Dita cy ee 
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where, 
Cn = =(. Ty p sinnpdp 
(9) 
4. SOLUTION OF THE PROBLEM 


The solution of the problem to be 
obtained is 
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Where, all conditions of convergence are same 
as in (2)[10]. 


Proof: 
Choose 


f(p) = u(p, 0) = To 


- oO- m, . (4)1),, 
=p 1(sin p) 1G ecsinp)?t | (11) 


1,q 


Now combining (11) and (9) and 
making the use of the integral (4), we derive 


Ch = 

w,\ (1 

> meme [CSE SaMoua, 
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(12) 


Putting the value of Cn from (12) in (8), 
we get the required result (10). 


5. SPECIAL CASES 
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On specializing, the parameters G-function may 
be reduced to several other higher 
transcendental functions [11]. Therefore, the 
result (10) is of general nature and may reduce 
to be in different forms, which will be apply in 
the literature on Applied Mathematics and other 
branches. 


CONCLUSION 


Finally, observed result [10] will be linear and 
have homogeneous boundary-conditions and 
the separation of variables method and Graph 
Theory can be applied to several types of 
problems on _ bounded _ three-dimensional 
domains. This result of this classical problem 
will be applied in three dimensions, the cooling 
of a sphere. The assumed symmetries in this 
problem will permit us to reduce the dimension 
of the problem to spatial dimension and time. It 
has also a great scope of future and new 
scholars. 
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